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Abstract 

The reactions pn — > duj and pn — > d<j) are studied within a relativistic meson- 
exchange model of hadronic interactions. Predictions for the total cross sec- 
tions and for the angular distributions of the vector mesons are presented. 
The resulting cross sections near threshold are around 10 - 30 fih for pn — > duj 
and 200 - 250 nb for pn — > dcj). A moderate deviation of the cross section ratio 
(Tpn^d<t>/ f^pn-^du) from that of the Okubo-Zweig-Iizuka rule is predicted. 

PACS: 13.75.-n, 14.20.Dh, 25.10+s, 25.40-h 
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I. INTRODUCTION 



The possibility that there might be a rather significant strangeness content in the nu- 
cleon, as indicated by the analysis of the TT-nucleon E term [1] and the lepton deep-inelastic 
scattering data by the EMC and successor experiments [2] , has triggered a wealth of exper- 
imental and theoretical investigations on this topic over the last decade. In this context, 
hadronic reactions at low and medium energy play an important role. Here one expects that 
the presence of strangeness in the nucleons would manifest itself in counting rates that sig- 
nificantly exceed predictions based on the Okubo-Zweig-Iizuka (OZI) rule [3,4]. Indeed, the 
experiments on antiproton-proton (pp) annihilation at rest at the LEAR facility at CERN 
revealed a strong violation of the OZI rule for various decay channels involving the meson 
[5-7]. Substantial deviations from the OZI rule were also found in the reactions pd -^^He0 
[8] and pp — > (jxp [9]. More recently the DISTO collaboration reported on an experiment 
on exclusive production of and u mesons in proton-proton (pp) collisions [10]. In this 
measurement a cross section ratio of about a factor 10 larger than the estimate from the 
OZI rule was observed (at a cms energy of about 80 MeV above the production threshold). 
Further experimental investigations of the reactions pp pp(f) and pp ppu, as well as 
pn —>■ d(p and pn dur, closer to their thresholds, are planned for the near future at the 
COSY facility in Jiilich [11,12]. 

Recently we have developed a model for the reactions pp — > ppcu and pp — > pp(f) within the 
conventional meson-exchange picture [13,14]. In this model meson production is described by 
the two mechanisms depicted in Fig. 1, namely the nucleonic current and the npv {v = uj, 0) 
meson exchange current. Those two processes were identified as providing the dominant 
contributions to these vector-meson production reactions. We employed this model for a 
combined analysis of close-to-threshold and cu production data with special emphasis on 
the apparent OZI rule violation reported by the DISTO collaboration. It turned out that 
the data do not require a large NN(f) coupling constant. Indeed, the range of values for gNN4, 
extracted from this analysis was found to be compatible with the OZI value [14] - which 
might be interpreted as an indication that there is no need for introducing an ss component 
into the nucleon in order to describe those data. On the other hand, it was necessary to 
introduce a violation of the OZI rule at the 7rp0 vertex in the meson-exchange current - 
which, in turn, provided the enhancement of the production cross section over the OZI 
estimation as seen in the DISTO data. 

In the present paper we apply this model to the reactions pn — > du and pn d(j). First 
of all, it is desirable to have quantitative results for those reaction channels. In this context 
we would like to emphasize that all our model parameters have ah^eady been fixed in the 
previous study [14], where they were constrained by requiring a consistent description of the 
reactions pp ppu and pp ppcj). Therefore we are able to provide genuine predictions 
for observables for pn — > du) and pn — > dep. Such predictions should be very useful for 
corresponding measurements that are planned at the COSY facility in Jiilich. Furthermore 
they can be used for testing the reliability of our model when data become available in the 
future. 

Secondly, in our previous studies [13,14] we found that the angular distribution of the 
produced vector mesons is very sensitive to the production mechanisms. We showed that 
this angular distribution enables us to determine both the relative importance of the basic 
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reaction mechanisms (cf. Fig. 1) and their absolute magnitTidcs. In fact, it was this special 
feature that enabled us to extract the (l)NN coupling constant from the analysis in Ref. 
[14] . Thus it is interesting to see whether a similar signature is also present in the reactions 
pn — > dcu and pn — > dcj) and, specifically, whether these reactions are even better suited for 
discriminating between the different production mechanisms. 

Finally, in the aftermath of the LEAR experiments, Ellis et al. [15] proposed a qualitative 
model in which it is assumed that the nucleon contains, besides u and d quarks, an intrinsic 
ss component. In this model 4> production proceeds via a "shake-out" or a "rearrangement" 
of the s and s quarks in the nucleon. Assuming, furthermore, that the ss component in 
the nucleon is polarized introduces selection rules for the various production channels and 
allows a qualitative interpretation of the state dependence shown by the pp annihilation 
data. Recently this picture has also been used for some qualitative estimates for the 
production in NN collisions [16]. It is interesting to explore whether the predictions of our 
conventional model differ from the results that would follow from the picture of Ellis et al. 

The paper is structured in the following way: In Section II we give a review of our 
relativistic meson-exchange model for vector-meson production. In Section III we present 
and discuss our results for the reactions pn — > dcu and pn d(j). The paper concludes with 
a short summary. 

II. THE MODEL 

We study the reaction pn dv {v — a;, 0) in a DWBA approach [13,14]. The transition 
amplitude is given by 

M^=<^/>,|J'^|0,> , (1) 

where ipd is the deuteron wave function and 0j is the four- component unperturbed NN wave 
function in the initial state. We take the deuteron wave function of the model Bonn B as 
defined in Table A.l of Ref. [17]. (Note that this model has been also used to generate the 
pp final state interaction in our previous study [13,14] of the pp — > ppv reaction.) J'* is the 
current describing the vector-meson production process. In our model is given by the 
sum of the nucleonic and vp7r meson-exchange currents, = J^^^ + J/^jg^, as illustrated 
diagrammatically in Fig. 1. In principle there are, of course, many more meson-exchange 
currents that could contribute. Their relevance was thoroughly investigated in [14] and it 
was found that all other exchange currents yield contributions that are much smaller than 
those of the vpn meson-exchange current. 

Explicitly, the nucleonic current is defined as 

JLc- E {r>;is,u + uis,r^) , (2) 

with denoting the NNv vertex and Sj the nucleon (Feynman) propagator for nucleon j. 

The summation runs over the two interacting nucleons, 1 and 2. U stands for the meson- 
exchange A^A^ potential. It is, in principle, identical to the driving potential [17] used in the 
calculation of the deuteron wave function ipd in Eq. (1), except that here meson retardation 
effects are retained following the Feynman prescription. Eq. (2) is illustrated in Fig. la. 
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The structure of the NNv vertex, (the subscript j = 1,2 is omitted), required in 
Eq. (2) for the production is obtained from the Lagrangian density 



where ^{x) and V^{x) stand for the nucleon and vector-meson fields, respectively. Qnnv 
denotes the vector coupling constant and Ky = Jmnv/ Gnnv-i with Jmnv the tensor coupling 
constant, tjin denotes the nucleon mass. 

As in most meson-exchange models of interactions, each hadronic vertex is furnished 
with a form factor in order to account for, among other things, the composite nature of 
the hadrons involved. In this spirit the NNv vertex obtained from the above Lagrangian is 
multiplied by a form factor. The theoretical understanding of this form factor is beyond the 
scope of the present paper; we assume it to be of the form 

^^^^'^ = A^^ + (Z^ - mlY ' 

where denotes the four-momentum squared of either the incoming or outgoing off-shell 
nucleon. It is normalized to unity when the nucleon is on its mass-shell, — nv^. 

The Vfm vertex required for constructing the meson-exchange current, J^g^ (Fig. lb), is 
derived from the Lagrangian density 

C,^{x) = -^^Sapu.d^f^ix) ■ d^T^{x)V^{x) , (5) 

where £ai3vii denotes the Levi-Civita antisymmetric tensor with £0123 = — 1- The vfm vertex 
obtained from the above Lagrangian is multiplied by a form factor which is taken to be of 
the form 

It is normalized to unity at = and = xm^. The parameter a; (= or 1) is introduced 
in order to allow for different normalization points in order to take into account that the 
values of the vpir coupling constants that we employ are extracted at different kinematics 
[14]. gr^p^ is determined at — rri^ and q^ — m^, whereas gu,pn is extracted at g'^ = and 
q^ — m^. Accordingly, we use the normalization x — 1 for the form factor at the former 
vertex and x = for the latter. 

The meson-exchange current is then given by 

nn-k{(1-k)]i + {1^2) , (7) 

where D^piqp) and A(g^) stand for the p- and 7r-meson (Feynman) propagators, respectively. 
The vertices F involved are self-explanatory. Both the NNp and A^A^tt vertices, F^^^ and 
Fatattt, are taken consistently with the A^A^ potential used to generate the deuteron wave 
function in Eq. (1). 
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Our model described above contains five parameters to be fixed for a given vector-meson 
produced: two for the mesonic current (the couphng constant g^p^^ and the cutoff parameter 
Am) and three for the nucleonic current (the coupling constants Qnnv and Jnnv, and the 
cutoff parameter Atv). 

The coupling constant g^^pTr was extracted from the measured decay width (0 — > p + tt) 
and ^i^pTT from the radiative decay width {uj ^ 'j + tt) assuming vector meson dominance. All 
other parameters were fixed in a combined analysis of all available near-threshold data on 
4>- and a;-meson production in pp coUisons [14] . Besides using the data from both reactions, 
we also assumed relations between corresponding parameters in the production amplitudes. 
Specifically, we assumed that the form factors at the (ppir and cupir vertices are the same. 
Likewise, we assumed that the form factor at the meson production vertices in the nucleonic 
current are the same. The NNu coupling constant is fixed to the SU(3) value, QNNiu = 9, 
based on Qnnp of Ref. [18]. Furthermore we assumed that also suggested 

by SU(3) symmetry, with —0.5 < k < 0.5. 

Details of our strategy for fixing the model parameters are outhned in Ref. [14]. Here 
we only want to mention that the lack of a more complete set of data prevented us from 
achieving a unique solution. Rather, we arrived at four different model solutions that all 
provide a comparably good description of the data but differ in the value of the nucleon 
cutoff mass Ajv, the value of 5fjv7v</> and the value of k, cf. Table I. We shall employ all these 
model solutions in the present investigation. 

Finally, we want to mention that the initial state interaction (ISI) is not taken into 
account explicitly in our model. Rather, its effect is accounted for effectively via an appro- 
priate adjustment of the (phenomenological) form factors at the hadronic vertices, cf. the 
corresponding statements in Ref. [14]. Admittedly, this procedure leads to a certain degree 
of uncertainty of our predictions because the ISI in the reactions pn — > dv and pp — > ppv, 
respectively, is governed by different partial waves (different isospin states). On the other 
hand, we don't have at present any better and more solid alternative for dealing with the 
ISI interaction. There arc no NN interaction models in the literature that arc still valid 
in the energy range relevant for uj and (p production for an explicit treatment of ISI effects. 
Furthermore, the qualitative prescription proposed in Ref. [19], based on phase shifts and 
inelasticity parameters, can't be used either because there is no reliable phase shift analysis 
beyond 1.3 GeV for the isospin T — channel [20] that we need here. 

III. RESULTS AND DISCUSSION 

The data that were used in Ref. [14] to fix the free parameters of our model consisted 
of some near-threshold total cross sections for the reaction pp — > ppui [21] and of the total 
cross section, as well as the angular distribution of the at Tiab — 2.85 GeV for the reaction 
pp ppcj) [10]. Naturally, the more selective experimental information available for the 
latter reaction (angular distributions) put also a stronger constraint on the model parameters 
and, as a consequence, the predictions for pn dcp for the four parameter sets in Table I 
don't differ much. As a matter of fact, the general features of the results are practically 
identical and therefore we show only the predictions for one of the sets (set 1 of Table I) in 
Fig. 2. Values for the total cross section at a few excess energies Q are, however, compiled in 
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Tabic II for all parameter sets. It is evident from Fig. 2 that the reaction pn dd is strongly 
dominated by the meson exchange current. The relative magnitude of the contribution from 
the nucleonic current is even smaller than in case of pp — > pp4>; it was about 18 % there, here 
it is a mere 3 %. The stronger suppression of the nucleonic current in the reaction pn — > d0 is 
due to the fact that the meson exchange current is weighted by a larger isospin coefficient (-3 
for pn, — > d(f) as compared to +1 in the case pp — > pp(l)). The nucleonic current, however, has 
contributions from the exchange of iso-scalar as well as iso- vector mesons. Consequently, only 
the latter are multiplied by a larger isospin coefficient. Furthermore, since the isospin factor 
changes sign, there are now also cancelations between the meson exchanges that contribute 
to the nucleonic current. As in the case of pp pp(t>-i there is a destructive interference 
between the nucleonic current and the meson exchange current which, however, is now much 
less pronounced. Nonetheless, the predicted angular distributions of the emitted meson, 
which are fairly flat, show some downward bending at very forward and backward angles 
as a consequence of this destructive interference between the two production mechanisms, 
cf. Fig. 3. For higher energies (Fig. 3, Q = lOOMey) the suppression of the angular 
distribution at forward and backward angles becomes more pronounced and it should be 
possible to verify experimentally whether the angular distribution does indeed have such a 
structure. In general there is not much difference between the angular distributions resulting 
from the four parameter sets, except that those with «; = 0.5 yield a total cross section which 
is about 20 % larger, cf. also Table II. 

Predictions of total cross sections for the reaction pn — > duj for the four parameter sets 
are shown in Fig. 4. Also here we observe that the meson exchange current is the dominant 
production mechanism. Its relative importance is likewise increased compared to pp ppoo. 
Indeed, now the meson exchange current dominates for all four parameter sets, whereas in 
the reaction pp — > ppui it is the nucleonic current that provides the dominant contribution 
for the parameter sets 2 and 4. Still, in all cases, the nucleonic current is large enough to 
introduce sizable interference effects. Consequently, we obtain a much stronger variation 
of the predicted total cross sections for pn du than for pn, dcf). The results vary by 
a factor of roughly 3, as can be seen from Fig. 4 and also from Table II. Note that there 
are particularly strong interference effects in case of parameter set 2 (Fig. 4b), where the 
resulting total cross section is significantly smaller than the individual contributions. As 
we shall discuss later, the stronger variations in the uj production cross section have also 
consequences for the cross section ratio R(j,/ui = o'pn^dip/o'pn^duj relevant for the comparison 
with the corresponding ratio that follows from the OZI rule. 

We observe that for Q < 20 MeV, effects from the finite width of the u meson become 
important. These effects are included in our calculation by folding the calculated cross 
section with the Breit-Wigner mass distribution of the u meson and can be seen in Fig. 4. 

Predictions for the angular distribution of the uj meson are displayed in Fig. 5. Naturally, 
here we see also larger variations between the predictions for the different parameter sets. 
But, since the meson exchange current is the dominant production mechanism for all sets, all 
the angular distributions show qualitatively a similar structure, i.e. they are rather flat at the 
lower energy {Q = 30 MeV) and exhibit a more-or-less strong downward bending at forward 
and backward angles at the higher energy {Q = 100 MeV). We would like to emphasize, 
however, that this behavior differs from the one in the reaction pp — > ppuo. In the latter 
reaction a clear enhancement of the angular distribution in forward and backward direction 
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was seen in the DISTO experiment [10] (at Tiab — 2.85 MeV), and such an enhancement for 
the reaction pp ppuo is also predicted by our model (parameter sets 2 and 4) - even for 
energies closer to the uj production threshold. 

As a general feature, let us also mention that the cross sections for the reactions pn — > dv 
are significantly larger than the ones for pp — > ppv near threshold. Apart from the isospin 
factors involved, such a behavior can be understood in terms of phase space arguments (the 
two body phase space opens up more rapidly close to threshold) and is well-known from 
other meson production reactions (see, e.g., Refs. [22,23]). 

We now turn our attention to the cross section ratio R^/^J = o-pn^d<j}/<^pn^dw for ^ com- 
parison with the prediction that follows from the naive application of the OZI rule [4] . The 
adjective "naive" refers to the fact that the cross section ratio is simply equated to the 
square of the ratio of the relevant coupling constants at the (f) and uj production vertices. 
The OZI rule plus SU(3) imply that these coupling constants are proportional to sina^ and 
cosa^, respectively, where stands for the deviation of the mixing angle from its ideal 
mixing value of 9(^ideai) — 35.3°. Thus the OZI estimate for the cross section ratio is simply 
given by -R^/c^ = tan^ [4]. With = 3.7° as extracted from the quadratic mass formula, 
this yields a value of R(j)/i^ = 4.2 x 10"'^. 

Our results for R^/^^ are presented in Fig. 6. The cross sections are evaluated at the 
same excess energy in order to minimize effects from differences in the phase space. We see 
that the largest values are predicted for excess energies close to threshold. E.g., at Q 10 
MeV we get values between 32 x 10~^ (set 2) and 9 x 10~^ (set 3). Thus, compared to the 
ratio that follows from the naive OZI rule, -R^/o; = 4.2 x lO^"', our models yield a moderate 
deviation of at most a factor of 8. At higher energies, Q = 80 ~ 100 MeV, the enhancement 
over the OZI estimate decreases to a factor of around 3 or so. 

So, what is the origin of this enhancement over the OZI estimate as predicted by our 
model? As already mentioned in the Introduction, our model involves an explicit violation 
of the OZI rule at the 7rp0 vertex in the meson exchange current. As discussed in detail in 
Ref. [14], it had to be introduced in order to achieve a simultaneous and consistent description 
within our model of the available data on the reactions pp ppu and pp — > particularly 
of the energy dependence of the total cu production cross section and the angular distribution 
of the meson. This explicit OZI violation, in terms of the 7rp0 and npu coupling constants 
used, suggests an enhancement of around 3 in the cross section ratio. With regard to the 
nucleonic current, the employed NNu and NN(f) coupling constants lead to results that 
exceed the OZI value only in one case, namely for parameter set 3, cf. Table I. The 
corresponding enhancement factor for the cross section ratio amounts to around 2. 

It is evident from Fig. 6 that the cross section ratios resulting from the full model cal- 
culation differ significantly from those values implied by the employed coupling constants. 
Obviously dynamical effects such as interferences, energy dependence of the production am- 
plitudes, etc. play a rather important role here. For example, parameter set 3 - which is the 
only model where the NNcf) coupling exceeds the OZI value - yields the smallest enhance- 
ment over the naive OZI estimate among all models considered, cf. the dash-dotted curve 
in Fig. 6. Furthermore, the particularly large enhancement over the OZI estimate found for 
parameter set 2 (especially at small excess energies) is, in fact, due to a suppression of the 
u! production cross section (resulting from strong interference effects as pointed out above) 
and not caused by an enhancement in the production, cf. Table II. Apparently fairly 
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large deviations from the naive OZI prediction can be generated by dynamical effects within 
the conventional picture: i.e., without introducing any "exotic" mechanisms. Consequently, 
one should be very cautious in drawing direct conclusions on the strangeness content in 
the nucleon from such cross section ratios. The behavior of the cross section ratio as the 
excess energy approaches zero is due to the finite width of the a;-meson which prevents the 
a;-meson production cross section from decreasing rapidly as it does in the case of 0-meson, 
cf. Figs. 2,4. 

Evidently, the enhancement of a factor of 3 or so over the OZI estimate at higher excess 
energies in Fig. 6 is much smaller than the enhancement over the OZI estimate of about a 
factor 10 found for o^^pp^japp^-p^ by the DISTO collaboration [10] at an excess energy 
Q ~ 80 MeV of the 0. However it is important to realize that their measurement was done at 
a fixed incident beam energy {Tiah = 2.85 GeV); therefore the corresponding excess energy of 
the produced lu is already 319 MeV. Although corrections for the differences in the available 
phase space were obviously applied when extracting the above result, there are other effects 
that influence the ratio such as the energy dependence of the production amplitude, the onset 
of higher partial waves, etc., which one cannot correct for easily. Therefore it is possible 
that the actual deviation from the naive OZI rule in the reactions pp ppv is also smaller. 
Indeed our results for (Jpp^pp^/ app^pp^^, shown in Fig. 7, support this conjecture. We see here 
that the results are very similar to the ones for crpn^d<j}/<^pn^dw Specifically, at the excess 
energy of Q = 80 MeV the enhancement over the OZI estimate is also only around a factor 
3. Thus it would be interesting to perform a measurement of the ratio (7pp_>pp<^/(Tpp_»ppa; at 
the same (or at least similar) excess energies. 

Finally, let us comment on a recent work by Ellis et al. [16]. This work builds upon the 
assumption that the nucleon wave function contains an admixture of negatively polarized ss 
quark pairs and allows one to deduce some qualitative predictions for the state dependence of 
(j) production in NN, but also in A^A^ collisions too. Specifically, this picture has been used to 
explain the observed differences in the branching ratios for pp — * (/)7r° at rest - which is rather 
large and OZI violating when the initial pp system is in a "^Si state, but small and consistent 
with the OZI rule for the ^Pi initial state. With regard to (p production in A"A^ collisions, 
it is argued that production should be likewise much larger from initial spin-triplet states 
[16]. Specifically, this means that the production cross section in the pp induced reaction 
is expected to be enhanced over the OZI rule because, close to threshold, it will proceed 
predominantly via the transition ^Pi — SqS. (We use here the standard nomenclature for 
labeling the partial waves of the A^A^ states and the angular orbital momentum of the 
relative to the final A^A^ system.) On the other hand, for pn induced reactions like pn — > 
the production amplitude is given by the transition ^Pi — >^ SiS and should therefore be OZI 
suppressed. In contrast, our dynamical model predicts cross section ratios of comparable 
magnitude for the pp — > ppv and pn dv reactions, as can be seen from a comparison of 
Figs. 6 and 7. Clearly, experimental information on these ratios would be very useful. 

IV. SUMMARY 

We have investigated the reactions pn — > d0 and pn — > du near threshold in a relativistic 
meson-exchange model. The elementary vector-meson production mechanism is provided by 
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the nucleonic and the 0p7r (and cupyr, respectively) mcson-cxchange currents. The presented 
results for total cross sections and for the angular distributions of the vector mesons are true 
predictions. All model parameters have been taken from a previous study of the reactions 
pp ppcj) and pp ppu [14] . 

A clear consequence of the model is that the reaction pn — > d(f) is almost completely 
dominated by the meson exchange current. The contributions from the nucleonic current, 
which were found to be already small for pp pp(p, are now even less significant - especially 
for the total cross section. The angular distribution of the 4> meson, however, still reflects 
some influence of the nucleonic current through interference effects. It is predicted to be 
rather flat for excess energies Q ^ 30 MeV but to show an increasing suppression of very 
forward and backward angles for higher energies. 

In case of the reaction pn duo the contributions from the nucleonic current are some- 
what larger. Thus, stronger interference effects between this current and the meson-exchange 
current are possible. As a consequence, the predictions for pn — > duj based on the four pa- 
rameter sets of Ref. [14] differ much more than those for pn — > d^, showing variations by 
up to a factor 3 in the total cross sections. This also implies greater variations in the cross 
section ratio (Tpn^d<f>l (^pn^dw, for which our model yields values that exceed the estimate 
based on the naive OZI rule by factors of up to 8. 
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TABLES 



TABLE I. Parameters of the four different models employed in the present study. The val- 
ues g^ppT^ = -1.64, g^pT^ = 10, gNNoj = 9, and Am = 1450 MeV are used for all four sets. 

K = Kv = fNNv/gNNv with V = (f),LLl. 

Set 1 2 3 4 



K -0.5 -0.5 +0.5 +0.5 

Ajv 1170 1411 1312 1545 

gNN4, -0.45 -0.19 -0.90 -0.40 



TABLE II. Total cross sections for the reactions pn — dw and pn — d(f) resulting from the 
model parameter sets 1 - 4 at selected excess energies Q. 
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FIGURES 



FIG. 1. (j) and w-meson production currents, J^, included in the present study: (a) nucleonic 
current, (b) vfm meson exchange current. M = tt, 77, p, w, cr, Oq. 

FIG. 2. Total cross section for the reaction pn d(j) as predicted by the model set 1 of Table I 
as a function of excess energy Q. The solid line is the result of the full model. The dashed 
(dash-dotted) curves show the contributions of the nucleonic current (meson exchange current) 
alone. 

FIG. 3. Angular distribution in the center-of-mass system of the produced (/)-meson in the 

reaction pn d(j) at the excess energies Q = 30 MeV (left panel) and 100 MeV (right panel). The 
dashed, solid, dashed-dotted and dotted curves are the results for the model parameter sets 1, 2, 
3, and 4. 

FIG. 4. Same as Fig. 2 for the reaction pn — > duj as predicted by the model sets 1 (a), 2 (b), 
3 (c), and 4 (d) of Table I. The dotted curves correspond to the results with the sharp mass of 
728.6MeF for the w-meson. 

FIG. 5. Same as Fig. 3 for the reaction pn doo. 

FIG. 6. Ratio of the total cross sections CFpn^d(t>/^pn^du as a function of excess energy Q. The 
dashed, solid, dashed-dotted and dotted curves are the results for the model parameter sets 1, 2, 
3, and 4. The horizontal solid line indicates the value predicted by the naive OZI rule. 

FIG. 7. Same as Fig. 6 for the ratio CTpp-^ppcf,/ CTpp-^ppuj- 
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